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In all known superconductors the pairing of fermions is not sensitive to the sign of their time- 
argument difference, leading to a Cooper pair wavefunction that is even in time/frequency. Four 
decades ago it was suggested that odd-frequency superconductivity should in principle be realizable. 
However, observation of odd-frequency superconductivity in bulk materials has remained elusive. 
Solving the field-dependent anisotropic Eliashberg equations, we present ab initio evidence for the 
emergence of odd-frequency pairing under an applied magnetic field in the archetypal superconductor 
MgE> 2 . We provide the full momentum, frequency and spin resolved dependence and magnetic field- 
temperature phase diagrams of the even and odd-frequency superconducting pair amplitudes and 
predict fingerprints of the odd-frequency state in tunnelling experiments. 


The pair amplitude of the two fermions forming the 
Cooper pair has to obey the Pauli exclusion principle, 
which implies that it has to be odd under particle ex¬ 
change. It is commonly assumed that the pair amplitude 
is invariant, that is, even, with respect to the sign of 
the fermionic time-argument difference, which leads to 
the standard classification of Cooper pair wavefunctions 
in terms of their symmetry under spatial and spin rota¬ 
tion [1]. It was however pointed out by Berezinskii [2] 
that the anti-symmetry of the pair amplitude can be ful¬ 
filled equally well when it is odd in time/frequency. Af¬ 
ter odd-frequency SC, which cannot be described within 
BCS theory, was first proposed for 3 He [2] it has since 
then been invoked in several occasions, including disor¬ 
dered Fermi liquids [3], High-T c compounds [4], heavy 
fermions [5] and hydrates [6]. Odd-frequency SC was also 
predicted to occur at superconductor to ferromagnet or 
to normal metal interfaces where the spatial symmetry is 
broken [7-10]. However, in spite of active research and 
the important ramifications for our understanding of su¬ 
perconductivity unequivocal observation of these exotic 
states remains elusive [11]. 

Remarkably, it was not until recently that the stability 
of bulk odd-frequency states was theoretically established 
[12, 13], thus opening up the possibility of experimental 
observations. Among all potential states, odd-frequency 
s-wave spin triplet (OST) SC [2] appears as the most 
robust since it can endure pair breaking by static non¬ 
magnetic impurities [4]. However, in a bulk material, 
OST SC cannot prevail over the even-frequency counter¬ 
part unless an extremely strongly coupled and retarded 
electron-phonon (EP) interaction is at play [14]. Un¬ 
til a material that fulfills these conditions is found, a 
promising alternative route is to look for special cases 
in which odd-frequency pairing could be possible [15- 
17]. The most plausible situation arises when an even 
frequency s-wave spin singlet (ESS) superconductor is 
placed in an external magnetic field, since the breaking 
of time-reversal symmetry is then expected to induce an 
OST SC component [16]. 

Despite recent progress, a fully microscopic, material- 


specific theory that supports the reality of OST SC is 
missing. We have extended the anisotropic Eliashberg 
framework to include the effect of an external magnetic 
field, allowing for OST SC as a selfconsistent solution, 
and show here the unambiguous existence of OST SC in 
the two-band superconductor MgB 2 - 

The two-band superconductor MgB 2 (T c = 39 K) [18, 
19] is an ideal candidate to study the magnetic field in¬ 
duced OST SC. MgB 2 is described thoroughly by ab ini¬ 
tio methods [20]. Its Fermi surface consists of sheets with 
primarily it or a orbital character [21]. High frequency in¬ 
plane boron modes lead to an enhanced, very anisotropic 
EP coupling that mediates the pairing [22, 23]. Near 
the center of the Brillouin zone (BZ) where these modes 
soften, the EP coupling well exceeds A q « 2 [23] (see also 
supplementary). Thus, for some regions in momentum 
space, the EP interaction in MgB 2 becomes both strongly 
retarded and coupled, satisfying the conditions for odd- 
frequency pairing. Due to the inherent anisotropy, SC 
in MgB 2 is characterized by a very anisotropic s-wave, 
two-gap structure (A^ = 2.8 meV and A a = 7meV) [24] 
which can be explained with unprecedented precision by 
fully anisotropic Eliashberg theory [25, 26]. Our ab ini¬ 
tio results for zero field completely confirm the previous 
calculations, which emphasizes the accuracy of our self- 
consistent solutions for finite magnetic fields. 

In Fig. 1 we show the calculated magnetic field- 
temperature phase diagram for the ESS and OST super¬ 
conducting components. For H > 0, the OST SC appears 
and coexists with the ESS SC in the same H-T regime 
where the latter is non-zero. The H-T dependence of the 
ESS SC follows that of a Pauli limited superconductor; 
it goes to zero monotonously via a second order phase 
transition with temperature, except for a very narrow 
region near the upper critical field (H p ), where the tran¬ 
sition becomes first order. The transition to the normal 
state with the field changes from first to second order 
above T^26K. Our solutions provide the first ab initio 
prediction for the paramagnetic limiting field of MgB 2 
(see, e.g. [27]) when orbital effects are neglected. We 
find H P =119T, which is significantly less than previous 
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FIG. 1. Temperature and magnetic field dependence of superconductivity in MgB 2 . a, H-T phase diagram of the 
even-frequency superconducting component. Dashed (solid) lines denote second (first) order transitions, b, H-T phase diagram 
for the odd-frequency superconductivity. The insets show the Matsubara frequency dependence of each component for several 
magnetic field values. The two gap structure, characteristic of SC in MgB 2 , can also be discerned, c, Band-resolved magnetic 
field dependence of the even and odd SC pairing amplitudes at low temperature. The lines correspond to the maximum values 
in Matsubara space of the momentum averaged superconducting pairing fields on each band, i.e. the peaks of the insets in a, 
b. Note that the OST pairing amplitude is an order of magnitude smaller than the respective ESS one. 


estimations (140 T) [28]. At low temperatures, the two 
ESS SC gaps remain almost constant up to ~ 39 T, 
as shown in Fig. lc. This is the lowest field strength 
needed to overcome the binding energy of the Cooper 
pairs originating from the 7r-band. Increasing the field 
further, leads to partial destruction of the latter gap and 
the concomitant decrease of the a gap, due to the strong 
interband EP coupling in this material. 

The OST component behaves similar to the ESS as 
concerns the order of the transitions to the normal state 
and its temperature evolution. However, it exhibits a dis¬ 
tinctively different magnetic field dependence; it follows 
a re-entrant behaviour with the field, as is shown in Fig. 
lb. The trend becomes clearer at higher temperature 
where the re-entrance peak moves away from H p . From 
Fig. lc one can observe that there exists, up to H « 39 T, 
a linear dependence of OST SC on the magnetic field ir¬ 
respective of the band index. Notably, the OST SC also 
has a two-gap structure but the pairing is much weaker 
than the respective ESS. 

Next, we examine the Fermi surface momentum de¬ 
pendence of the SC gaps, presented in Fig. 2. Since 
in the Eliashberg framework the superconducting gap 
retains its full frequency dependence (cf. Figs. 2C,F), 
it is more meaningful to define a superconducting gap 
edge, instead. In the zero-field case this is given by: 
ReAk(^k) = co>k [25]. However, here the Zeeman effect 
lifts the degeneracy between spin a =t, i quasiparticles. 
As a consequence, external probes see a superconducting 
gap edge that is effectively split in two. Moreover, the 
presence of the OST SC further modifies the gap edge. 
The relevant expressions can be found from the poles in 


the system’s Green’s functions at the Fermi level: 
w k,<r = —<7H(k, w k>(T ) + A e (k, w k(T ) + crA 0 (k, w kcr ) (1) 

where A k '°) is the order parameter for ESS (OST) SC 
and H k , w is the renormalized magnetic field term that 
includes selfenergy contributions (see supplementary). In 
the above, all quantities are real and the OST spins are 
perpendicular to H. We define an effective gap edge for 
each SC component as: A^^(k) = A e ( 0 )(k, cj k ^))- 
Since the real part of the OST order parameter van¬ 
ishes at cj = 0 (Fig. 2f), the gap edge of a pure odd- 
frequency superconductor should be zero. However, due 
to the presence of the ESS SC and the magnetic field, 
this is finite in our case. In Fig. 2a we show the full mo¬ 
mentum dependence of the spin-^ ESS gap edge over the 
Fermi surface of MgB 2 . The 3D tubular networks are due 
to the 7r-bands, while the almost two dimensional cylin¬ 
ders are due to the cr-bands [21]. The two-gap structure 
is similar to the zero-field case [25] but the momentum 
anisotropy differs. For the 7r-bands we find gap values 
between 1.8 — 3.5 meV, which are close to the zero-field 
values, while for the cr-bands we find gap values between 
5.6 — 8 meV and thus an enhanced anisotropy. Fig. 2b 
shows the difference between the spin-split gap edges of 
the ESS component. Remarkably, the splitting is more 
efficient for the cr-bands where the gap edge is larger. 

As shown in Fig. 2d, the OST SC gap edge is very 
anisotropic and at H = 40 T it is two orders of magnitude 
smaller than the respective ESS. Comparing Figs. 2D and 
2B we find that AJ(k) is proportional to the difference 
of the ESS gap edges for the two spin components. This 
finding directly evidences that the OST superfluid den¬ 
sity is proportional to the number density of spin flipped 
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FIG. 2. Momentum and frequency dependence of superconductivity in MgEU at T = 4.2K and H = 40T. a, 

The Fermi surface of MgB 2 in the conventional Brillouin zone coloured by the values of the even-frequency superconducting 
gap edge for spin-| quasiparticles. The distribution of the gap edge values is shown in the inset below, b, The difference 
between spin-^ and spin-t even-frequency superconducting gap edges, (c,) Typical frequency dependence of the real (Ag(w)) 
and imaginary (A'Jfa)) part of the ESS pairing amplitude at a given point on the Fermi surface, d, e, Same as in a, b but 
for the odd-frequency superconducting gap edge, f, Same as in c but for the real (Ao(u;)) and imaginary (A^w)) parts of the 
OST pairing amplitude. Near uj = 0, A ' 0 {uj) increases linearly with frequency whereas A"(c j) is almost constant. The latter is 
finite at uj — 0 as shown in the inset. 


carriers that participate in the ESS Cooper pairs. There¬ 
fore, the OST is subordinate to the ESS SC. 

From the spectral function of our system, we have 
derived the quasiparticle DOS in the superconducting 
state, N(uj) = \^ (7 N ct {uj). This quantity is propor¬ 
tional to the differential conductance measured in tun¬ 
nelling experiments [26]. The spin-resolved quasiparticle 
DOS reads 

= Re /_ |^ + o-H k ,J _\ (2) 

N f \ [ (w + - (A£ >u + aA° > J2] i 

where Np is the DOS at the Fermi level in the normal 
state and (.. .)k denotes Fermi surface averaging. Us¬ 
ing our selfconsistent data in Eq. (2) we calculate the 
magnetic field evolution of the tunnelling spectra at low 
temperature, shown in Fig. 3a. For H = 0 we find two 


peaks around 3 meV and 7 meV that signal the quasipar¬ 
ticle excitations above the i r and the a gaps, respectively 
[25, 26]. We also obtain the recently observed fine mo¬ 
mentum structure of the latter gap [26]. As the field 
increases, the peaks begin to split due to the spin degen¬ 
eracy lifting. This is clearly visible already at 10 T for the 
7r peak, while at 20 T the splitting in both peaks should 
be very clearly resolved. The gap starts to close around 
39 T where the magnetic field strength is comparable to 
the minimum gap edge value on the 7r-bands. 

The odd-frequency component is finite for finite field 
and increases linearly with uj in the low frequency regime 
(cf. Fig. 2f). Therefore, one would generally expect the 
OST gap to manifest as a shift of the tunnelling peaks for 
finite cj, to higher (lower) uo due to destructive (construc¬ 
tive) interference with the dominant ESS gap. However, 
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FIG. 3. Dependence of single-particle tunnelling spectra on the magnetic field, a, Low temperature superconducting 
density of states of MgB 2 for several values of the external magnetic field, b, Comparison between tunnelling spectra when the 
odd-frequency pairing is included in the Eliashberg calculation and when it is not (black lines). As the field increases, the odd- 
frequency gap induces a shift in the tunnelling peaks of the order of 0.1 meV. This shift is pronounced for the superconducting 
gap over the cr-bands. c, In the absence of external bias, we predict a non-zero signal in the tunnelling spectra that is 
proportional to the imaginary part of the OST SC. This has a distinct magnetic field dependence which detection will serve as 
definite proof for the existence of odd-frequency SC in MgE> 2 . 


due to the smallness of the OST gap, these effects are 
practically washed out by the dominant ESS component 
except for two regions. The first is for fields near 39 T, 
just before the i r gap begins to close. There, the OST 
gap reduces the magnetic field threshold of the zero en¬ 
ergy peak by 0.5 T. The second is at frequencies near the 
cr-band peak, where the OST component shifts the peak 
to the right. For 40 T the calculated peak shift is around 
100 fieV and increases with the field as seen in Fig. 3b. 

We now focus on the zero bias regime, u) = 0. This case 
is of special interest, since all lifetime effects coming from 
the EP and the Zeeman interaction are absent. The real 
part of the induced OST gap is also zero (Fig. 2f) but the 
situation differs from that in a pure OST superconductor 
since here the ESS component provides a robust gap at 
the Fermi level. Hence, no low lying excitations should 
be expected at low temperatures and small, comparing 
to the 7r gap, magnetic fields. This is certainly true for 
H = 0 (Fig. 2a). However, when the field is turned on, 
the imaginary part of the OST component, A"(k, 0), be¬ 
comes finite (Fig. 2f). This term effectively contributes 
an additional lifetime broadening to the quasiparticles 
and remarkably enforces a nonzero result in Eq. (2) even 


at cj = 0. If the OST is not included in the theory, 
the latter is ideally zero. Therefore, a zero-bias signal 
in the single-particle tunnelling appears due to the in¬ 
duced OST superconducting component. This signal has 
a distinct magnetic field dependence, as is shown in Fig. 
3c. This zero-bias plateaux is even more pronounced for 
T<4.2K where the OST SC is more favored (Fig. lb). 
Thus, we propose that a measurement of the zero-bias 
tunnelling under magnetic fields less than ?^39T and at 
low temperatures would serve as definite proof for the 
identification of odd-frequency SC in MgB 2 and in any 
superconductor that supports this kind of state. 


Our precise quantitative results predict the existence 
of OST SC in MgB 2 , provide an in-depth insight into 
the microscopic nature of this exotic state and pave 
the way for the ultimate experimental identification of 
odd-frequency SC. Furthermore, our approach initiates a 
novel ab initio “roadmap” for the search of such exotic 
magnetic-field-induced phenomena relevant for both bulk 
and interface physics [7, 9, 10]. 
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